Dynamical conductivity of ungated suspended graphene 
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Frequency dependent conductivity of Coulomb interacting massless Dirac fermions coupled to 
random scalar and random vector potentials is found as a function of frequency in the regime 
controlled by a line of fixed points. Such model provides a low energy description of a weakly 
rippled suspended graphene. The main finding is that at the neutrality point the a.c. conductivity 
is not frequency independent and may either increase or decrease with decreasing u, depending on 
the values of the disorder variances A^, Aa and the Coulomb coupling a = e 2 /(evp). The low 
frequency behavior is characterized by the values of two dimensionless parameters 7 = A^/a 2 and 
Aa which are RG invariants, and for small values of which the electron- hole "puddles" are effectively 
screened making the results asymptotically exact. 



The physics of massless Dirac fermions in two spatial 
dimensions has received renewed attention since their dis- 
covery in single-layer graphene 0, 2,3]. The great interest 
is not unrelated to the quantum critical^, IE Hj S B Hi 
nature of the system near the neutrality point, where the 
Fermi level lies precisely at the (Dirac-like) band crossing. 
Indeed, absence of any intrinsic long distance lengthscale 
sets constraints on the (low) frequency or temperature 
dependence of any physical quantity. In this regard, elec- 
trical conductivity a plays a special role since in two spa- 
tial dimensions it is expected to be proportional to e 2 /h; 
the proportionality constant, which need not be finite, 
depends only on the nature of the renormalization group 
(RG) fixed point characterizing the low energy-long dis- 
tance physics 0, [E EE EL EH- Electrical conductivity 
measurements at the neutrality point therefore consti- 
tute a direct probe of the non-trivial physics emerging at 
the end of the RG trajectory. 

Recent experiments performed on monolayer graphene, 



14, 



If 



have 



both suspended [13J and on the substrate 
found that the optical conductivity near the neutrality 
point is <t(uj) = J 4-, i.e. largely frequency indepen- 
dent and equal to ir/2 it the natural units. What lit- 
tle frequency dependence there is in the regime where 
Huj ~ leV can be attributed to the curvature correc- 
tions to the electronic dispersion which deviates from 
the perfectly conical massless Dirac-like at such large 
energies [IE]. While presently there is no conductivity 
data at the lower frequencies of interest (at sub meV 
scales) , it is natural to ask whether such frequency inde- 
pendent cr(w) should persist down to huj ~ fc^T. Since 
the role of charged impurities located at the substrate 
is naturally eliminated in the suspended samples, the 
dominant source of scattering is most likely the random 
configuration of strain due to the graphene sheet rip- 
pling and possibly from the boundary effects imposed by 
the scaffolding necessary for the actual suspension. Such 
long wavelength strain fields are known to couple to the 
massless Dirac particles of grap hene as a vector poten- 
tial and a scalar potential IE 12, IE], and their combined 



effect on the a.c. conductivity, together with the effects 
of the electron-electron (Coulomb) interactions, are ana- 
lyzed below. 

In the non-interacting model of massless Dirac par- 
ticles coupled to the random scalar and random vector 
potentials with variances and A a, it has been long 
known that within the perturbative RG, A^ grows upon 
approaching low energies and the theor y fl ows to a per- 



turbatively inaccessible fixed point [13, [20] . The effects 
of Coulomb interactions, parameterized by a dimension- 
less coupling a — e 2 /evp, and random vector potential 
(without scalar potential randomness) on the a.c. con- 
ductivity has been studied in Ref.[3] where it was found 
that the a.c. conductivity is non-universal and depen- 
dent only on A a which is marginal in the RG sense. Such 
non-universality is directly tied to the appearance of the 
infra-red (IR) locally stable line of fixed points. The two 
loop effects on the RG flow diagram have been incorpo- 
rated in Ref . [2l| and extended to other types of disorder 
in Ref [H]. 

In this work, we study the combined effects of the 
unscreened Coulomb interactions, the quenched random 
vector and scalar potential disorder which arise naturally 
in the model of (randomly) strained sample, and ana- 
lyze the frequency dependent conductivity within per- 
turbative RG. This includes the effects of monolayer rip- 
ples and electron-hole "puddles". The perturbative RG 
adopted here is technically much simpler than the large 
N approximation adopted by Foster and Aleiner 22j to 
map out the phase diagram, and to leading order, leads 
to qualitatively similar results (for differences beyond 
the leading order, and the advantages of the former, see 
Ref.[2lj]). Moreover, the weak coupling RG can be easily 
extended to the calculation of cr(w), which was not cal- 
culated in (22^ . The main finding is that cr(u>) is not fre- 
quency independent and may either increase or decrease 
with decreasing u>, depending on the values of A^, A a 
and a. The low frequency behavior is characterized by 
the values of the dimensionless parameters 7 — A^/a 2 
and Aa which are RG invariants. The ultimate low fre- 
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FIG. 1: The renormalization group flow diagram in the scalar 
disorder A^ - Coulomb interaction a = e 2 /(evF) plane. 
There are two marginal parameters (RG invariants): the 
variance of the random vector potential Aa and the ratio 
7 = A^/a 2 . For 7 < 32/(7t 2 Aa) each RG trajectory can cross 
three fixed points (the strong coupling one is not shown) . The 
middle one, which is IR stable, merges with the IR unstable 
at a multicritical point (red circle) when 7 = 32/(7t 2 Aa). 
ao = 4Aa I tx lies along the fixed line discussed in [J, |21( . For 
7 > 32/(tt 2 Aa), there are runaway flows with no perturba- 
tively accessible fixed points. The phase diagram splits natu- 
rally into two regimes: Regime I (shaded) given by the locus 
of points which eventually run into the perturbatively accessi- 
ble IR stable fixed line show above, and Regime II (unshaded) 
with runaway flows. The u — > limit of the (collisionless) a.c. 
conductivity along the IR stable fixed line is given in Eq.JTJ, 
and its u dependence is discussed in the text. 



quency behavior depends on being in one of two regimes 
(see FigEJ. 

Specifically, for the bare couplings in the Regime I, 
which is determined by the conditions a < 8A A /^ and 
7 < tt 2 /(32A a ) or a > 8A A /n and A < (tt/2)q - 2A, 
then as oj — > 



e 

Ti 



6 



(23 - 67r)7r 2 
96j 



32 A 



We start with the imaginary time partition function 

Z = J v^e' (Sa+Sd " +SM) (2) 

where 

So = [ dr I d 2 ri/)(r, t) (9 T + vf<J ■ p) ^>(r, t) (3) 

Sdis = 



Si, 







dr I d 2 rip(r, r) (0(r) + v F a ■ a) tp(r, r) (4) 

11 

1 f 13 

dr ld 2 rd 2 r'^(r, T )V(\r - r'\)iptp(r' ,t)(5) 





The last term corresponds to the (Coulomb) electron- 
electron interaction V(\r — r'\) — £ | r 1 r / | , where e is the 
dielectric constant which may differ from 1. We assume 
that the disorder is uncorrelated with variances: 



4k4') = (2^) 2 5(k-k> 2 A 
i(X,} = (2TT) 2 5(k-k')S^A A . 



(6) 
(7) 



As has been discussed extensively in the past, the scalar 
and vector potentials are naturally connected to the 
appearance of strain tensor ug as 4> = g(u xx + U yy ), 
a x = b(u yy - u xx ), a v = 2bu xv H, Oj, EE [23Lwith the 
estimates g « 20 - 3CW0, [jj] and b w A^pl, [23 1. 

We can perform the (quenched) average over the gaus- 
sian disorder fields <fi(r) and a(r) using the standard 
replica trick of including n copies of the fermion fields: 
ip — * ip'\ where i = 1, 2, . . . , n. The resulting replica field 
theory is 



(Z n ) di , 



Z)i(j i i(; i e~( So + S '>' + SA+Si ' lt ' ) 



(8) 



subject to the constraint < 7 < ir 2 / (32Aa). In this pa- 
rameter regime the weak coupling RG flow equations lead 
to a perturbatively accessible IR stable line of fixed points, 
making the above result asymptotically exact. The ac 
conductivity, in the collisionless limit of interest here, 
therefore behaves as a universal amplitude, depending 
only on the RG invariants A^ and 7. In the Regime II, 
which encompasses the parameter regime not included 
in Regime I there are no perturbatively accessible fixed 
points, and the problem remains open. The details of 
the calculation leading to the above claims, as well as ui 
dependence of the conductivity, are presented below. 



where 
So = 

= 

(i) s A = 

Sint = 



dr / d 2 ^ (d T + v F a ■ p) i> 1 



v 2 F A q 



drdr' J d 2 v^ l ^\r,TW^(r,T') 
drdr'f d 2 r^V M i/; i (r, t)^ (r, r') 



2 Jo 



dr / d 2 rd 2 rVV l (r, r)V(\r - r'|)ff(r',T) 



As usual, we assume large momentum cutoff A for the 
above fermion modes, and perform the renormalization 
of the bare coupling constants e 2 ,vf -i^^/^^- To first 
order we find that the imaginary time Greens function 
(Fig.©) satisfies 



V F 



(9) 
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FIG. 2: The self energy diagrams to leading non-trivial order 
in {A^, Aa} (left) and a. (right). 



The renormalization condition demands that we absorb 
the dependence on the cut-off A into a field rescaling con- 
stant Z and the bare couplings. We do so at an arbitrary 
scale uj = k = k where we demand that 



Giu,(k)|„ = ZG£(k)|„ = Z(-iK + v$a ■ k)- 1 



(10) 



The renormalized Greens function G^(k), at any to and 
k, is now independent of A. This leads to the RG equa- 
tion for the Fermi velocity 



Pvf — 



dvi 



d log A 



A A A e 2 \ 



2tt 



4v F ) 



(11) 



FIG. 3: Diagrams which contribute to the disorder renormal- 
ization. 

To determine the RG scaling of the disorder vari- 
ances we need to analyze the p functions of the ef- 
fective replica coupling constants. We do so by writ- 
ing the equations for the irreducible four point vertex: 
rW = g| = (r^^^)con.amp. The renormaliza- 
tion prescription demands that at a scale uj = k = k, 
r(4) = ^ T f. This means that at arbitrary a;,k, the 
quantity Z 2 T^ can be made independent of A. To this 
order in coupling constants (Fig. ([3])) we find 



Pa, = 
Pa a = 

Pe* = 



dA 
d log A 

dA A 
d log A 

de 2 



A0 
2tt 



A. 



4vf 



= 0, 



(12) 
(13) 
(14) 



d log A 

which agrees with Ref . ■ As argued in Ref . 0] the last 
equation is exact. The corresponding flow diagram is 
shown in Fig[TJ 

Defining the dimcnsionlcss Coulomb coupling constant 
a = e 2 /evp, the above equations imply the existence of 
two RG invariants 7 = A^/a 2 and A a, i.e. 



dA, 



d log A d log A 



0. 



(15) 



FIG. 4: Diagrammatic contribution to the (a.c.) electrical 
conductivity within Kubo formula. From left to right: free 
Dirac fermion contribution, disorder vertex and self energy 
corrections, Coulomb interaction vertex and self energy con- 
tribution. 



Since conductivity does not acquire anomalous dimen- 
sion we have 



d 



dlogA 

where the differential operator 

d „ d 



B ) a(uj; A, A^, A^, v F , e ) = 0, (16) 



B = Pa, 



dA d 



Pa, 



dA, 



P. 



d 

V F "a 
OVf 



Pe 



de 2 ' 



(17) 



The solution of the above RG equation must satisfy the 
scaling law 



a{uj- pA, A 4> {pA),A A {pA),v F {pA),e 2 {pA)) 
a(uj; A, A4A),A A (A),v F (A),e 2 (A)) 



(18) 



where p is a positive real number [24|. 

A pedestrian perturbation theory calculation to the 
leading order in coupling constants (FigHJ gives 



cr P t(uj) = 4- 



A^ 
24 



Aa 
24 



10 



25 
~6~ 



(19) 



This extends the result found in [4[ to include the scalar- 
disorder potential contribution. Note that the above ex- 
pression does not satisfy the scaling law (fT8|) . since ac- 
cording to Eqs. (|lHll4|) the coupling constant a and the 
scalar disorder variance A^ do have non-trivial depen- 
dence on A. Nevertheless, to the same order in the cou- 
pling constants, both (fT8|) and (fT9|) can be satisfied if 



a(uj) = 4- 



24 



Aa 
24 



16 



25 



(20) 



where we used A c f ) (p) = ja 2 (p). The dimensionless 
Coulomb coupling constant a{p) is defined as the solution 
of 



da 
dlogp 



, 1 2 
' 1 to" 



Aa 

7T 



with the initial condition ot(l) = 
pendence can be found implicitly 



(21) 

a. The functional de- 



« 2 (f) 



+ A) ( 

\ A 



TV 
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%a-(l-A) VO-1) 



It, 



A 2 



-a— A(A —1) 

! (22) 
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FIG. 5: Illustrative dependence of a on frequency for Aa = 1 

2 

and 7 = 0.95 t^-. The lower curve corresponds to the start- 
ing condition a m (near the repulsive fixed point at the 
origin of Fig.fTJ) and the upper to ot as 3.25 (near the strong 
coupling repulsive fixed point). Notice that since both curves 
asymptote to the same fixed point, the u> — > limit of a is 
the same and the value is given by Eq.{T]). 



where A = wl — ^■yA A . The above equation is easily 
inverted numerically. As an illustration, the resulting 
conductivity as a function of frequency for A A = 1 and 
in the vicinity of the multicritical trajectory is plotted in 
Fig ©. 

The explicit dependence on u> can be found in some 
limiting cases. In the vicinity of the IR stable fixed line, 
but away from the multicritical point (see Fig[T]), we find 



47 i _ 



^" 47a \ V ^ 7T 2 7^ 



where the (crossover) exponent 



I67 



1 



32 



7A A 1 



32 



jA A 



(23) 



(24) 



Note also that 
d log uj 



4 



7T 
16 



25 



7 /k\ 

12 a( A } 



da 
d log uj ' 



The Eq.(j2Tl) then implies that for initial a < j 1 - - n) , 
the conductivity in the Region I to the left of the IR fixed 
line (light grey shaded portion of Region I in FigJTJ in- 
creases with decreasing uj. On the other hand, to the 
right of the IR fixed line (dark shaded portion of Re- 
gion I) a(uj) decreases with decreasing uj. In cither case, 
however, as uj 0, a asymptotes to the value given by 
Eq.([T]). Moreover, while in the latter case, the large fre- 
quency limit is outside of the scope of perturbative RG, 



in the former case the high frequency limit of conductiv- 

Thus, the appearance of the infra-red (locally) stable 
line of fixed points at finite Coulomb coupling and finite 
disorder (Fig. [T]) provides a natural theoretical avenue to- 
wards non-universality of the longitudinal electrical con- 
ductivity, since the precise position along such line is typ- 
ically beyond experimental control. Nevertheless, any 
correlation between the variance of the (independently 
measurable) strain configurations and the minimal con- 
ductivity would provide a good test of the above theory. 
It is also important to address the effects of the general 
form of disorder. As argued in Ref. 22j, if one starts with 
the most general disorder potential allowed by the sym- 
metry of the graphene honeycomb lattice, the RG flow 
trajectories runaway to strong coupling/strong disorder 
and the physics is perturbatively untractable. The pic- 
ture presented here may nevertheless be physically rel- 
evant for suspended samples, since the primary mode 
of coupling of the smooth deformations to the graphene 
Dirac fermions is via scalar and vector potentials, and all 
but the absence of the short wavelength components of 
the scattering potential gives the above (globally unsta- 
ble) fixed line a large basin of influence. 
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